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LONG  TERM  GOALS 

The  principal  long-term  goals  of  this  project  are  to  provide: 

•  In  concert  with  related  projects,  an  improved  understanding  of  vertical  turbulent  mixing 
and  it  consequences  in  stratified  estuarine  flows. 

•  Improved  methods  for  determination  of  the  suspended  particulate  matter  (SPM)  concentra¬ 
tion  and  transport  fields  tfom  acoustic  and  optical  backscatter  measurements. 

OBJECTIVES 

The  past  year’s  work  focused  on  the  second  goal,  finding  improved  methods  for  determination  of  the 
SPM  concentration  and  transport  fields  tfom  acoustic  and  optical  backscatter  measurements,  given  ap¬ 
propriate  calibration  data.  As  a  part  of  this  effort,  we  have  used  continuous  wavelet  transform  methods 
to  determine  the  importance  of  tidal  and  non-tidal  mechanisms  in  horizontal  SPM  transport. 

APPROACH 

In  coordination  with  the  parent  grant  of  this  AASERT,  we  have  developed  single  and  multiple- 
tfequency  inverse  analysis  methods  for  SPM  determination.  The  student  supported  by  this  AASERT 
for  1999-2000,  Annika  M.  V.  Fain,  focused  on  the  single  tfequency  inverse  techniques.  A  single- 
tfequency  inverse  analysis  uses  dynamical  information  present  in  SPM  profiles.  Profiles  (“basis  func¬ 
tions”)  are  defined  for  each  time  each  of  a  small  number  of  size  or  JVg  classes.  The  contribution  of  each 
basis  function  to  the  observed  profile  at  a  given  time  is  determined  by  an  inverse  analysis  (via  non- 
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negative  least-squares  or  NNLS).  Results  depend,  however,  on  the  assumed  SPM  dynamie  balance 
used  to  form  the  basis  functions.  Previous  analyses  have  used  a  balance  in  the  vertical  between  turbu¬ 
lent  SPM  flux  and  particle  settling.  Fain  (2000)  and  Fain  and  Jay  (2000)  also  use  a  modified  Rouse  bal¬ 
ance  approach,  but  have  improved  upon  previous  single-sensor  analyses  in  several  respects.  Fbst,  the 
presence  of  aggregates  of  unknown  density  means  that  scattering  behavior  and  settling  velocity  (Ws) 
are  not  simple,  known  functions  of  particle  size.  When  aggregates  are  present,  we  formulated  the 
analysis  around  settling  velocity  (Ws)  classes  rather  than  size  classes.  Second,  the  non-negative  least 
squares  (NNLS)  inverse  technique  used  is  quite  sensitive  to  mismatch  between  the  assumed  and  actual 
settling  velocity  distribution.  We  have,  therefore,  used  known  Ws  spectra.  Finally,  our  inverse  tech¬ 
nique  is  based  on  a  broader  dynamical  analysis  that  defines  estuarine  SPM  dynamics  in  terms  of  five 
non-dimensional  numbers.  The  parameters  are  the  Rouse  number  P,  an  advection  number  A  (a  non- 
dimensional  height  of  maximum  SPM  concentration  above  the  bed),  a  trapping  efficiency  E  (measuring 
the  efficiency  of  retention  of  material  in  the  system),  a  trapping  potential  Tp  (ratio  of  stratification- 
induced  shear  to  shear  velocity),  and  a  Supply  Number  Sr  {Sr  is  a  measure  of  the  strength  of  the  fluvial 
or  marine  SPM  supply  to  the  system).  These  parameters  have  been  chosen  both  for  dynamical  rele¬ 
vance,  and  because  they  are  relatively  robust  against  calibration  difficulties. 

WORK  COMPLETED 

The  single-frequency  code  has  been  developed,  tested  and  used  in  analyses  of  the  SPM  dynamics  of  the 
Columbia  River  estuary.  This  work  forms  the  basis  of  Annika  Fain’s  M.S.  thesis  and  related  publica¬ 
tions  (Fain,  2000,  Fain  et  al,  2000,  Jay  et  al,  2000).  The  single-frequency  code  is  also  used  in  the 
multi-frequency  inverse  analysis  described  in  the  report  for  the  parent  grant  of  this  AASERT. 

RESULTS 

Analysis  of  four  7-8  mo.  ADCP  records  (from  the  ONR-supported  CORIE  system)  has  allowed  for  ex¬ 
ploration  of  the  above  parameter  space  (Eain,  2000;  Eain  et  al,  2000,  Jay  et  al,  2000).  This  data  set  is 
favorable  for  a  test  of  the  inverse  method,  because  a  very  large  calibration  data  set  is  available,  and  in¬ 
terference  from  zooplankton  in  the  acoustic  signal  was  small.  An  important  result  is  the  dependence  of 
trapping  efficiency  E  on  river  flow  (Eigure  1).  E  in  the  Columbia  River  estuary  is  shown  as  a  function 
of  the  seasonal  trajectory  of  Supply  Number  Sr  for  the  Columbia  River  for  May  to  December  1997. 
This  period  encompassed  the  strongest  freshet  since  1974,  and  1997  had  the  highest  total  annual  flow 
of  any  year  in  the  20^*'  century.  Also  shown  is  E  for  the  Eraser  River  estuary  under  extreme  freshet  con¬ 
ditions  in  July  1999,  during  one  of  the  strongest  freshets  of  the  20^*'  century.  Results  suggest  that  the 
moderate  flows  after  the  spring  freshet  resulted  in  the  most  efficient  trapping  in  the  Columbia.  E  was 
very  low  in  the  Eraser,  because  of  the  high  flow.  These  results  conform  to  intuitive  understanding.  If 
the  river  flow  is  too  high,  trapping  of  SPM  in  the  estuary  will  be  small  (low  E)  despite  a  very  high  sup¬ 
ply,  because  the  estuary  is  so  short  that  all  salt  (and  SPM  retained  in  the  salt  wedge  salinity  intrusion) 
will  be  removed  on  each  greater  ebb.  If  the  flow  is  too  low,  the  estuary  is  very  large,  but  upstream  bot¬ 
tom  currents  are  weak  and  may  not  retain  SPM  effectively  (low  E  again).  Also,  the  limited  SPM  supply 
is  diluted  in  a  very  large  estuarine  water  mass.  E  is  then  maximized  at  intermediate  flow  levels  where 
upstream  bottom  flows  are  robust,  but  salt  is  not  removed  on  every  ebb.  It  appears  that  the  1997  data 
set  encompassed  the  period  of  maximum  E  for  the  Columbia.  The  residence  time  of  river  water  and 
SPM  in  the  Eraser  was  so  small  in  July  1999  that  SPM  was  exported  over  the  top  of  the  salt  wedge  be¬ 
fore  it  could  settle  into  the  lower  layer,  yielding  E  «\.  There  is  one  other  important  influence  not  con¬ 
sidered  in  the  two-dimensional  (x  and  z)  dynamics  encompassed  by  the  parameter  space  described  here 
-  SPM  can  be  stored  in  peripheral  areas  during  freshets  and  supplied  to  the  main  body  of  the  estuary 
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even  during  periods  of  low  river  flow.  This  faetor  likely  influenees  the  value  of  Sr  that  shows  maximal 
E  in  Figure  1 .  The  asymptotie  low  and  high  Sr  behaviors  posited  are,  however,  likely  robust. 

The  single-sensor  inverse  approaeh  deseribed  here  has  the  advantage  of  dynamieal  and  numerieal  sim- 
plieity,  but  there  are  two  primary  weaknesses:  a)  negleet  of  horizontal  adveetion  of  SPM,  and  b)  the 
faet  that  a  bulk  SPM  ealibration  is  employed,  whereas  baekseatter  strength  depends  on  partiele  size. 

The  method  is,  then,  best  suited  to  eireumstanees  where  the  range  of  sizes  present  is  not  too  large.  The 
presenee  of  aggregates  (the  dominant  material  found  in  the  Columbia  River  estuary)  is,  however,  not  an 
obstaele. 

IMPACTS/APPLICATIONS 

Applieation  of  this  teehnique  should  allow  the  routine  interpretation  of  aeoustie  baekseatter  (e.g.,  from 
aeoustie  Doppler  eurrent  profiles)  in  terms  of  SPM  properties,  provided  that:  a)  the  baekseatter  field  is 
dominated  by  SPM  not  zooplankton,  and  b)  adequate  ealibration  data  are  available. 

TRANSITIONS 

Two  eontinuous  wavelet  transform  tidal  analysis  programs  with  sample  data  sets  illustrating  their  use 
have  been  plaeed  on  the  Pi’s  web  page  (see  CWT  software  library  sub-heading  under  http://www.ese. 
ogi.edu/~djay).  Dr.  R.  Signell  of  the  US  Geologieal  Survey  has  also  set  up  a  link  to  these  programs  on 
his  Sea-Mat  web  site  (http://erusty.er.usgs.gov/sea-mat).  These  programs  and  related  eodes  have  been 
used  in  the  analysis  of  the  SPM  data  derived  from  inverse  analyses. 

RELATED  PROJECTS 

Work  for  the  Tidal  Channels  AASERT  projeet  has  been  eoordinated  with  the  National  Seienee  Foun¬ 
dation  Columbia  River  Land-Margin  Eeosystem  Researeh  (LMER)  Program  and  with  the  Oeeano- 
graphie  and  Environmental  Charaeterization  of  Coastal  Regions  (OECCR)  funded  by  ONR. 
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Figure  1:  The  May-December  seasonal  trajectory  in  the  Columbia  River  (CR)  estuary  of  SPM  parti¬ 
cle  trapping  (T rapping  Efficiency  E)  as  a  function  of  SPM  supply  (Supply  Number  Sr),  A  point  for 
the  Fraser  River  (FR)  is  added  to  demonstrate  the  high  river-flow  asymptote. 
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